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SECTION  1 
INTRODUCTION 


THE  PROBLEM  ADDRESSED  BY  CEP 

This  paper  describes  an  algorithm  for  calculating  Circular  Error 
Probability  (CEP).  CEP  relates  measurement  errors  to  the  errors  in  a 
calculation  that  is  computed  from  those  measurements.  As  such,  it  is 
a  useful  tool  for  assessing  the  accuracy  and  sensitivity  of  any  cal¬ 
culated  quantity  to  etrored  inputs.  CEP  is  frequently  used  with 
target  location  estimates.  It  is  obtained  under  the  assumption  that 
the  measurement  variances  are  known. 

THE  INTERPRETATION  OF  CEP 

The  CEP  associated  with  an  estimated  target  location  (x,y)  is 
defined  to  be  the  radius  of  the  smallest  circle  with  center  at  (x,y) 
which  has  a  50X  probability  of  containing  the  true  target  coordinates. 
Note  that  the  value  of  CEP  Is  a  distance  associated  with  the  50X 
probability  and  not  a  probability  itself,  as  is  suggested  by  the  name. 
CEP  provides  a  measure  of  the  search  area  within  which  the  true  target 
location  can  be  expected  to  be  found.  The  better  the  accuracy  of  the 
(x,y)  estimate,  the  smaller  the  value  of  CEP. 

ORGANIZATION  OF  THE  PAPER 

In  section  2  the  relevant  equations  for  CEP  computation  are 
developed.  Computational  steps  are  delineated  along  with  required 
assumptions  and  a  mathematical  foundation  is  provided.  A  simplified 
algorithm  that  approximates  CEP,  and  which  is  often  quoted  in  litera¬ 
ture,  is  discussed  in  this  paper.  The  algorithm  has  a  high  degree  of 
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accuracy.  A  full  derivation  of  the  algorithm  is  found  in  appendix  A. 
Also  included  in  section  2  is  a  discission  of  the  Trace  of  the  covar¬ 
iance  orror  matrix.  Under  certain  conditions  the  square  root  of  the 
Trace  is  shown  to  be  proportional  to  the  CEP. 

Section  3  is  a  summary  of  the  algorithmic  procedure  for 
obtaining  CEP  from  measurement  data.  It  is  the  vorking  guide  to  be 
used  in  CEP  computation.  Section  A  gives  a  brief  description  of  a 
planar  bistatic  environment.  This  geometry  is  frequently  used  in 
emitter  location.  Section  S  gives  examples  of  CEP  calculations  and 
their  interpretation  in  a  bistatic  environment.  Three  commonly 
occurring  cases  are  discussed.  Section  6  is  a  summary  of  the  paper. 


SECTION  2 

THE  HATHBHATICS  OP  CEP 


THE  FUNDAMENTAL  EQUATION 

By  definition  the  CEP  is  the  radius  of  a  circle  about  the  (x,y) 
estimate  which  has  a  SOX  probability  of  containing  the  target.  The 
equation  for  any  circle  of  radius  CEP  and  center  (x,y)  is 

(x'-x)2  +  (y'-y)2  -  CEP2  . 

Let  u  ■  x'-x  and  v  »  y'-y,  and  rewrite  the  circle  equation  as 

u2  +  v2  «  CEP2  . 

If  we  assume  (u,v)  are  Gaussian  variables  with  mean  (0,0)  and 
variances  ( au 2 , ov 2 ) ,  their  joint  distribution  is  given  by 


Hx)’-  H-r)1]  • 


and  the  requirement  that  the  CEP  has  a  50X  probability  of  containing 
the  target  is  expressed  by  the  integral  relation 


To  obtain  the  correct  value  of  CEP  one  must  solve  equation  (2). 
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Evaluation  of  (2)  involves  numerical  approximation  since  (except 
for  some  specific  values  of  ou  and  ov)  the  integral  has  no  known 
antiderivative.  Nevertheless,  an  explicit  formula  for  CEP  can  be 
obtained  vhich  approximates  the  solution  to  (2)  vith  a  high  degree  of 
accuracy.  An  approximating  formula  for  CEP  is  described  in  the 
following  subsection. 

A  SIMPLIFIED  FORMULA  FOR  CEP 

Whenever  the  following  two  conditions  hold, 

1.  the  error  variables  (dx,dy)  are  uncorrelated,  and 

2.  the  error  probability  densities  are  zero  mean  Gaussian 
distributions  with  variances 

then  the  following  approximating  formula  for  CEP  may  be  used: 


where 


CEP 


o,  <0.67+0.8v2) 

L 

0.59eL(l+w) 


if  0<  w  <0.5 

if  0.5<  w  <  1 


et  -  Maximum!  <JdJt,odyl  , 
es  »  Minimum! od)(,ody]  , 

w  «  W 


A  derivation  of  (3)  is  found  in  appendix  A. 


(3) 
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Quality  of  the  Approximation 


Figure  1  shows  the  graph  of  (CEP/oL)  versus  as/oL  for  both  the 
"exact"  value  of  CEP  and  the  approximation.  The  "exact"  value  was 
computed  from  equation  (2)  to  four  decimal  figure  accuracy.  The 
approximation  was  obtained  from  formula  (3).  A  measure  of  the 
difference  between  the  graphs  is  shown  by  figure  2. 

In  figure  2,  deviations  between  the  two  graphs  are  displayed  in 
terms  of  the  percent  error,  defined  by 


Percent  Error 


100 


<CEP>.x.ct  -  <CEP>. 


pp  r  ox 


<CEP>.x.ct 


Figure  2  shows  that  the  formula  (3)  approximation  differs  from  the 
true  CEP  by  a  maximum  of  1.5  percent,  with  the  average  err.»r  less 
than  one  percent.  The  approximation  is  thus  seen  to  be  an  extremely 
accurate  measure  of  true  CEP. 


DETERMINATION  OF  THE  STANDARD  DEVIATION 


Assumptions 

Use  of  the  CEP  formulas  require  knowledge  of  the  standard 
deviations  (<*dx>ffdy)-  1°  this  section  we  show  how  to  obtain 
<odx,ody)  in  terms  of  the  measurement  errors.  The  following 
assumptions  are  used: 

.  If  m  is  a  measurement  used  in  the  calculation  of  x  or  y, 
then  dm  (i.e.,  the  error  in  m)  is  assumed  to  be  a  Gaussian 
random  variable  with  mean  zero  and  with  known  standard 
deviation  a.  . 

a  m 
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CE PI  a 


Figure  2.  Percent  Error  in 
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.  The  standard  deviation  for  each  of  the  dm  variables  is 
sufficiently  "small"  so  that  second-order  terms  in  the 
Taylor  expansion  of  dm  can  be  neglected. 

.  All  of  the  dm  variables  are  statistically  independent. 

For  application  of  the  CEP  formulas,  one  must  be  sure  that  these 
assumptions  are  met. 

An  Expression  for  the  (dx,dy)  Errors 

Our  initial  objective  is  to  express  the  errors  (dx,dy)  in  terms 
of  the  dm  errors.  Initially,  ve  must  have  (x,y)  expressed  analyti¬ 
cally  in  terms  of  measurements.  If  this  can  be  done,  then  a  good 
approximation  to  the  (dx,dy)  errors  is  provided  by  the  derivative  of 
(x,y).  A  required  assumption  for  the  approximation  is  that  each  of 
the  dm  errors  are  "small."  For  illustration  ve  use  the  case  where  x 
and  y  are  computed  from  three  measurements  (m,  m,  m, )  as  we  trace 
through  the  mathematical  development.  The  extension  to  an  arbitrary 
number  of  measurements  follows  a  similar  pattern. 

The  quantities  (dx,dy)  are  estimated  by  taking  the  total 
derivative  of  x  and  y,  namely, 


*■  ■  §;*■! 

dy  .  d.t 


*  fe*. 

*  fe,*, 

(4a) 

*  I/** 

*  to/"*  • 

(4b) 

Relationship  (4)  shows  that  errors  dx  and  dy  are  expressed  as  a 
linear  combination  of  the  measurement  errors.  From  statistical 
theory,  the  linearity  provides  an  important  observation,  that  is  to 


8 


say,  if  each  of  the  error  variables  (dm1 ,dm2 ,dm3 )  are  assumed  to  be 
Gaussian  with  mean  (0,0,0),  then  the  variables  (dx,d'»)  will  also  be 
Gaussian  with  mean  (0,0).  Thus,  our  assumption  of  Gaussian  behavior 
for  the  measurement  errors  allows  us  to  describe  the  errors  (dx,dy) 
by  means  of  a  Gaussian  analysis. 

In  matrix  form,  the  (dx,dy)  relations  are  expressed  compactly  as 


'dx' 

dy. 


dmxl 

dm2 

Ldm3 


(5) 


where  M  is  the  Jacobian  of  the  transformation,  expressed  as 


‘  3x 

3 x 

3x 

3m3 

3m2 

3m3 

UL 

2l 

iL 

3m2 

3m3 

An  Expression  for  the  Covariance  of  (dx,dy) 


The  covariance  of  (dx,dy)  is  a  2x2  matrix  from  whose  elements 
the  standard  deviations  of  dx  and  dy  can  be  derived  (a  squired  next 
step  in  our  CEP  computation).  In  appendix  B  it  is  shown  that  the 
covariance  of  (dx,dy)  is  expressible  in  terms  of  the  matrix  M  and  the 
covariance  of  the  (dm2 ,dm2 ,dm3 )  variables  by 


Cov(dx,dy)  *  M  Cov(dml,dm2,dm3)  M*  ,  (7) 

where  H*  is  the  transpose  of  matrix  H.  Multiplication  of  the 
expression  in  (7)  gives 


Cov(dx,dy) 


aii 

2 

a21 

a22 

(8) 
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where 


ayv  * 

(f,)V 

*  (tjv 

(9a) 

aJi  “  (Im,) 

* 

(fj 

<9b> 

&JV  ♦ 

*  • 

(9c) 

Extracting  Standard  Deviations  from  the  Covariance  Matrix 

The  standard  deviations  (<xdx,ady)  can  be  extracted  from  the 
covariance  matrix  (8)  using  the  following  considerations. 

If  an  appropriate  coordinate  system  were  chosen  for  the  (dx,dy) 
joint  distribution,  and  (dx,dy)  were  independent  variables  (as  we  are 
assuming),  then  the  covariance  of  (dx,dy)  would  identify  with  the 
standard  deviations  (adx,ody)  by 


Cov(dx,dy) 


(10) 


In  general,  however,  the  a12  term  of  the  covariance  matrix  (8)  will 
not  be  zero  and  a  diagonalization  of  the  matrix  (equivalent  to  a 
rotation  of  the  axes)  is  required  to  convert  matrix  (8)  into  the  form 
of  equation  (10).  After  rotation,  equation  (8)  will  have  the  form 


Cov(dx,dy) 


0  ‘ 

\  . 


(ID 
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where  are  eigenvalues  of  (5).  The  eigenvalues  are  given  by 


Trace 

2 


(12a) 


x  Trace  \ 

Trace2 -4Det 

-  2  -  • 

2 

where: 

Trace  «  a 21  +  a22  and 

Det  -  a:ja22  a12 

• 

(12b) 


('2c) 


Since  eigenvalues  of  a  matrix  are  invariant  under  rotation,  the 
eigenvalues  in  (12)  must  equal  the  variances,  namely, 


<W  =  ("dx2,*dy2)  • 

As  a  caveat,  note  that  when  the  eigenvalues  are  computed  from 
(12),  it  is  not  clear  (without  further  analysis)  which  eigenvalue 
associates  with  o  2  and  which  associates  with  e.  2.  However,  this 

x  ay 

knowledge  is  not  required  for  computation  of  the  CEP.  Only  the  ratio 
of  smaller  to  larger  eigenvalues  is  needed  (see  CEP  formula  (3)). 


Significance  of  the  Trace  of  the  Covariance  Matrix 

Under  certain  conditions,  the  CEP  can  be  shown  to  be  proportional 
to  the  square  root  of  the  Trace.  The  behavior  of  the  Trace  thus 
provides  both  a  qualitative  and  quantitative  description  of  CEP 
behavior. 
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The  necessary  condition  involves  the  radicand  in  (12)  that  is 
used  to  obtain  the  eigenvalues  \  and  Xj.  For  simplicity,  ve  rewrite 
(12)  in  the  form 


XA  «  Trace  [  1  k.-  4Det_/Trace2  ^  (i-1(2)  .  (13) 

Vhile  the  positiveness  of  the  covariance  matrix  always  guarantees 
that  the  radicand  is  positive,  i.e., 

0  <  4Det/Trace2  <  1  , 

we  consider  the  cases  for  which  the  Det  is  significantly  smaller  than 
the  square  of  the  Trace,  i.e.,  the  condition 

0  <  4Det /Trace J  «  1  .  (14) 

For  this  condition,  we  can  neglect  the  second  term  in  the  radicand  of 
(13)  and  replace  it  by  zero  as  a  first  approximation.  The  two 
eigenvalues  (equivalent  to  the  two  variances  adx2  and  <Jdy2)  take  on 
the  value  X^Trace  and  Xj*0.  Thus,  the  larger  standard  deviation 

becomes  oL=^lTrace  and  the  smaller  is  ffs=0.  Substituting  these 
values  into  the  CEP  formula  (3),  we  obtain  w=0  and 


CEP  =  0.67  'Tfrace  .  (15) 

Therefore,  under  condition  (14)  the  square  root  of  the  Trace  can 
be  used  as  a  "quick  look"  estimate  of  the  CEP.  Some  examples  are 
discussed  in  section  5. 
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SECTION  3 

A  SUMMARY  OF  THE  COMPUTATIONAL  STEPS  INVOLVED 
IN  A  CEP  DETERMINATION  OF  TARGET  ACCURACY 


This  section  provides  a  working  guide  tor  computation  of  CEP. 

The  reader  is  reminded  that  underlying  assumptions  must  be  met  to 
insure  both  validity  and  accuracy  of  the  CEP  estimate.  Chief  among 
these  are  the  following: 

.  The  measurements  used  in  the  computation  of  target  location 
have  errors  that  can  be  treated  as  Gaussian  random  variables 
with  zero  mean  and  known  standard  deviation. 

.  The  measurement  errors  are  statistically  independent. 

.  The  standard  deviation  for  each  of  the  measurement  errors  is 

sufficiently  small  so  that  the  second  order  terms  in  the 
Taylor  expansion  of  the  meas'.rement  error  can  be  neglected. 

Under  these  conditions  the  computational  steps  for  obtaining  CEP  are 
summarized  below. 

Steps  in  the  Computation  of  CEP 

Step  1.  Establish  an  analytic  relationship  between  (x,y)  and 
the  measurement  variables. 

Step  2.  Obtain  partial  derivatives  of  (x,y)  with  respect  to  the 
independent  measurement  variables  being  used  in  the 
(x,y)  computation. 
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Step  3.  Substitute  the  partial  derivatives  of  (x,y)  into 
equation  (9)  to  obtain  the  matrix  components 

an*a2a*ai2* 

Step  4.  Substitute  the  matrix  components  into  formula  (12c)  to 
calculate  the  Trace  and  Determinant  (Det).  If 
relationship  4Det<<Trace7  holds,  then  the  product 

0.6?*’lTrace  provides  an  excellent  estimate  of  the 
CEP.  Alternatively,  continue  on  and  compute  the 
eigenvalues  \  and  Aj  from  formulas  (12a)  and  (12b). 

Step  5.  Obtain  the  standard  deviations  odx  and  edy  by  taking 
the  square  root  of  A1  and  Aj . 

Step  6.  Substitute  the  values  of  adK  and  ody  into  formula  (3) 
to  obtain  the  CEP. 


SECTION  4 

BISTATIC  TARGET  LOCATION 


In  bistatic  detection,  two  platforms  are  used  to  acquire  infor¬ 
mation  about  a  target  T.  The  platforms  may  be  either  moving  or 
stationary.  Either  one  (or  both)  may  be  illuminating  the  target.  In 
mathematical  analysis,  the  tvo  platforms  and  target  are  geometrically 
portrayed  as  a  triangle  on  which  simple  equations  can  be  used  to 
describe  distances  and  angles. 

GEOMETRY 

Figure  3  describes  a  planar  view  of  the  geometry  used  in 
bistatic  calculations.  Platform  A  is  located  at  the  origin  (0,0). 
Platform  B  is  located  at  (Rif0).  The  target  T  is  located  at  (x,y). 
All  these  coordinates  are  measured  relative  to  the  local  origin  A. 

In  practice,  the  line  directions  and  locations  of  the  A,B,T 
triangle  are  also  measured  relative  to  an  external  reference  frame. 
For  example,  if  platform  B  was  functioning  as  an  information  Collec¬ 
tor,  it  would  be  natural  to  collect  all  directional  information 
relative  to  its  own  heading.  In  this  paper  ve  assume  that  informa¬ 
tion  is  measured  relative  to  some  reference  direction  on  platform  B, 
and  that  this  direction  lies  in  the  same  plane  as  the  triangle  formed 
by  A,  B,  and  T. 

Figure  3  displays  the  angles  measured  relative  to  platform  B.  4t 
defines  the  bearing  angle  (measured  counterclockwise)  from  platform  B 
direction  to  the  target;  defines  the  angle  (measured  counter¬ 
clockwise)  that  the  line  connecting  the  two  platforms  makes  with  the 
platform  B  direction.  Similarly,  ©  -  defines  the  bearing  angle 
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from  A  to  the  target  T  as  seen  froa  tha  platfora  B  perspective.  8  is 
Measured  clockwise  froa  the  bearing  line.  Proa  the  geometry,  the 
base  angles  of  the  triangle  vith  vertices  A,B,T  are  given  by  0  +  8  - 
and  -  ♦t. 

TYPICAL  PLATFORM  MEASUREMENTS 

Various  combinations  of  inforaation  can  be  used  to  acquire  an 
estimate  of  target  location.  Some  of  the  common  measurements 
acquired  are  listed  belov. 

Rt  The  measured  distance  between  platform  B  and  platform  A. 

R1  has  an  associated  measureaent  error,  dRi,  whose  average 
value  is  zero  and  whose  standard  deviation,  ed|t  ,  is  knovn. 

i 

T.ie  measured  angle,  relative  to  a  platform  B  frame  of 
reference,  which  defines  the  direction  of  the  line  from 
platform  A  to  platform  B.  has  an  associated  measurement 
error,  d^,  whose  average  value  is  zero  and  vhose  standard 
deviation,  cd^  ,  is  knovn. 

♦t  The  measured  angle,  relative  to  a  platform  B  frame  of 
reference,  which  defines  the  bearing  of  the  line  from 
platform  B  to  the  target.  has  an  associated  measurement 
error,  d+t,  vhose  average  value  is  zero  and  vhose  standard 
deviation,  cd4  ,  is  knovn. 

t 


0  Th«  measured  angle,  relative  to  e  platfora  B  frame  of 
reference,  vhich  defines  the  bearing  of  the  line  from 
platform  A  to  the  target.  8  has  an  associated  measurement 
error,  d0,  whose  average  value  is  sero  and  vhose  standard 
deviation,  edt,  is  known. 

Rp  The  measured  distance  from  platform  A  to  the  target.  Rp 
has  an  associated  measurement  error,  dRp,  whose  average 
value  is  zero  and  vhose  standard  deviation,  edR  ,  is  known. 

p 

Rt  The  measured  distance  from  platform  B  to  the  target.  Rt 
has  an  associated  measurement  error,  dRt,  vhose  average 
value  is  zero  and  vhose  standard  deviation,  ed|)  ,  is  known. 

t 

AT  The  Time-Difference-of-Arrival  (TDOA)  of  a  signal  received 
by  platform  B  from  tvo  different  paths.  If,  for  example, 
platform  A  were  illuminating  a  target  (e.g.,  in  the  case  of 
a  radar),  or  interrogating  a  target  (e.g.,  in  the  case  of  a 
beacon  system),  the  signal  would  be  received  directly  by 
platform  B  (along  the  line  connecting  A  and  B)  in  time  T2 , 
and  received  again  at  a  later  time  T2  (along  the  lines 
connecting  A  to  T  and  T  to  B)  from  reflection  off  or 
response  by  the  target.  The  time  difference  would  be  AT  * 
Tj  -  Tx.  If  the  signal  is  a  response  from  the  target,  the 
time  T2  includes  any  fixed  and/or  random  time  delays 
introduced  by  the  target  transponder.  Note  that  in  the 
case  of  beacon  systems,  the  "target"  could  be  the  Interro¬ 
gator  and  A  the  transponder.  AT  has  an  associated  measure¬ 
ment  error,  dAT,  vhose  average  value  is  zero  and  vhose 
standard  deviation,  od6T,  is  known. 
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SECTION  5 
EXAMPLES 


Assuming  that  the  angle  ^  and  the  distance  Ri  betveen  the 
platforms  A  and  B  have  already  been  obtained  (i.e.,  the  base  of  the 
triangle  with  vertices  A,B,T  has  been  determined  in  the  geometry  of 
figure  3),  then  two  additional  measurements  relating  to  target  are 
required  to  compute  the  target  location.  We  describe  here  the  common 
cases  where  two  angle  bearings,  or  one  angle  bearing  with  TDOA 
measurement,  are  available. 

CASE  1.  TWO  BEARING  ANGLES  TO  TARGET  ARE  MEASURED 

Bearing  angles  9  and  are  acquired  with  known  standard 
deviations  ode  and  ,  respectively.  Geometrically,  target 

t 

location  (x,y)  will  be  the  intersection  of  the  relative  bearing  lines 
defined  by  II  +  0  -  ^  and  (see  figure  4).  Analytically, 

(x,y)  is  obtained  by  solving  the  simultaneous  equations 

Tan(9-^)  =  y/x  and  TanC^-4^)  =  -y/(Ri-x)  . 

The  partial  derivatives  of  (x,y)  with  respect  to  the  measurement 
variables  are  given  in  appendix  C.  The  Trace  and  Det  are  obtained  by 
substituting  these  derivatives  into  equation  (13).  Evaluation  gives 


Trace 


„  2  2  n  2  2 
Rp  ”de  *  Rt  °d*t 


n  2  2  .  2  ,  .  2 

*  Ri  'd*,  *  s,n  <*-*i)odRi 


sin2  (a) 
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Figure  4.  Case  1  Geometry:  Target  Location  by  the  Intersection 

of  Two  Bearing  Angles 
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and 


Det 


RpR2t<wV  »defld*1 


) 


Risin2 


( W<od0  *  <’d*1)0B1 


sin1  (a) 


where  a  «  ^  -  0. 


Figure  5  is  a  plot  of  CEP  contours  for  the  Case  1  geometry  with 
«  0  deg,  Ri  =  150  nmi  and  with  standard  deviations  ode  *  5  deg, 
ad+  =  5  deg,  ad^  *  0  deg,  odR  =  0  deg.  Platform  locations  A  and  B 

t  i  i 

are  on  the  axes.  The  graph  is  to  be  interpreted  as  follows:  Select 
a  point  (x,y).  Linearly  extrapolate  (as  a  first  approximation) 
between  curves  to  estimate  the  curve  passing  through  (x,y).  The 
contour  number  associated  with  that  curve  corresponds  to  the  CEP  for 
point  (x,y) . 


As  discussed  earlier,  the  square  root  of  the  Trace  may  be  used 
to  qualitatively  describe  the  CEP.  For  the  geometry  of  figure  4, 
division  by  zero  occurs  in  the  Trace  under  three  sets  of  conditions: 
1)  +t  .  0(  e  -  0,  2)  -  180,  e  =  180,  and  3)  -  90,  0  «  90. 

Conditions  1  and  2  occur  along  the  horizontal  axis  of  figure  4.  This 
axis,  and  a  small  band  around  it,  is  commonly  called  the  "black  hole" 
region.  Along  this  line,  the  CEP  is  infinite.  Targets  near  the 
horizontal  axis  are  located  very  poorly  using  this  measurement  set. 
Condition  3  corresponds  to  targets  being  extremely  far  from  platforms 
A  and  B.  Such  targets  will  also  have  extremely  high  CEP. 

In  general,  the  triangulation  process  of  Case  1  for  obtaining 
target  location  applies  to  any  two  bearing  measurement  systems  (e.g., 
two  ground  radars,  a  ground  radar  and  aircraft  radar,  etc.).  An 
example  of  Case  1  application  is  the  cooperative  triangulation  of  two 
E-3  aircraft. 
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CASE  2.  ONE  BEARING  ANGLE  TO  TARGET  AND  TDOA  IS  MEASURED  FROM 
ILLUMINATION  OF  PLATFORM  A  BY  THE  TARGET  (HYPERBOLA) 

The  geometry  for  this  case  is  described  in  figures  6a  and  6b. 
Bearing  angle  6  (or  +t)  is  acquired  with  known  standard  deviation  <xde 
(or  ad^  ).  The  time  difference  of  arrival,  AT,  of  the  target  signal 

t 

from  two  different  directions  is  also  acquired  at  point  B.  The  two 
directions  traversed  by  the  acquired  signals  are  as  follows: 

Direction  1  is  the  direct  distance  Rt  from  the  target  to  point  B. 
Direction  2  is  the  distance  traveled  by  the  signal  (e.g.,  radar  or 
beacon)  from  the  target  to  the  platform  A  plus  the  distance  from  the 
platform  A  to  platform  B  (i.e.,  R  +  RA).  Time  along  direction  1  is 
received  at  Tj .  Time  along  direction  2  is  received  at  T2 .  The  time 
difference  is  AT  =  T2  -  Tx .  If  the  signal  is  a  response  from  tne 
target,  the  time  T2  includes  any  fixed  and/or  random  time  delays 
introduced  by  the  target  transponder.  AT  has  an  associated 
measurement  error,  dAT,  whose  average  value  is  zero  and  whose 
standard  deviation,  ®d6T,  is  known. 

One  can  show  from  analytic  geometry  that  the  target  location 
will  lie  on  a  hyperbola  defined  by  AT  and  Rt  (see  figure  6b).  The 
hyperbola  has  foci  at  (-C0 AT/2,0)  and  (Ri+Co AT/2,0) ,  and  eccentricity 
equal  to  Ri/C)JAT,  where  Co  is  the  speed  of  light.  Geometrically, 
target  location  will  be  the  intersection  of  the  relative  bearing  line 
n  +  0  -  ^  (or  ^  -  +t)  with  the  hyperbola.  Analytically,  target 
location  (x,y)  is  obtained  by  solving  the  simultaneous  equations 

Tan(0-^)  »  y/x  (  or  Tan(^t-^)  =  -y/(Rt-x)  ) 

and 

R,  +  R.  -  R.  ■  C  AT. 

P  1  L  O 
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Bearing  line  vith  a  Hyperbola 


The  partial  derivatives  of  (x,y)  with  respect  to  the  measurement 
variables  are  given  in  appendix  C.  The  Trace  and  Det  are  obtained  by 
substituting  these  derivatives  into  equation  (13).  Evaluation  gives 


2Rpi(ode +ad^)(1_cos<x)  +  Co2ffdAT  *  [1+cos<*t  -♦i>l‘adR^ 


,2  _  2 


Trace  = 


(1-cosa)' 


and 


p<ffde  +  V^at  *  CUcos (♦t-»i»-  o^j 

(1-cosa)2 


2  2 


2  2 


Det  = 


where  a  =  $t  -  0. 

Figure  7  is  a  plot  of  CEP  contours  for  the  Case  2  geometry  with 
=  0  deg,  Ri  =  150  nmi  and  with  standard  deviations  cde  =  5  deg, 

=0  deg,  od4T  =  10-6  sec.  Platform  locations  A  and  B  are 

i 

marked  on  the  axes.  Similar  to  Case  1,  the  graph  is  to  be  inter- 
prt  sd  as  follows:  Select  a  point  (x,y).  Linearly  extrapolate  (as  a 
fi’-o.t  approximation)  between  curves  to  estimate  the  curve  passing 
through  (x,y).  The  contour  number  associated  with  that  curve 
corresponds  to  the  CEP  for  point  (x,y). 

As  discussed  earlier,  the  square  root  of  the  Trace  may  be  used 
to  qualitatively  describe  the  CEP.  In  the  above  equation,  division 
by  zero  occurs  when  cos(a)  =  1.  This  condition  occurs  along  the  hori¬ 
zontal  axis  of  figure  5  to  the  left  of  A  (when  0  =  180,  =  180), 

and  to  the  right  of  B  (when  0  =  0,  =  0).  For  this  geometry,  these 

external  intervals  on  the  horizontal  axis  (and  a  small  band  around 
them)  are  called  "black  hole"  regions.  Along  these  intervals,  the 
CEP  is  infinite.  Targets  near  these  intervals  have  extremely  poor 
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location  measurement.  In  contrast,  targets  between  A  and  B  along  the 
horizontal  axis  are  able  to  be  located  vith  good  accuracy. 

Tactical  Air  Navigation  System  (TACAN)  is  an  example  of  Case  2 
application.  Here  the  airborne  target  interrogates  the  ground  station 
(i.e.,  the  ground  station  at  location  A  is  the  transponder). 

CASE  3.  ONE  BEARING  ANGLE  TO  TARGET  AND  TDOA  IS  MEASURED  FROM 
ILLUMINATION  OF  THE  TARGET  BY  PLATFORM  A  (ELLIPSE). 

The  geometry  for  this  case  is  described  in  figures  8a  and  8b. 
Bearing  angle  9  (or  +t)  is  acquired  vith  known  standard  deviation  <xd6 
(or  od4  ).  'fhe  time  difference  of  arrival,  AT,  of  the  target  signal 

t 

from  two  different  directions  is  also  acquired  at  point  B  with  stan¬ 
dard  deviation  od4T.  However,  as  opposed  to  Case  2,  the  directions 
traversed  by  the  acquired  signals  are  different.  Direction  1  is  the 
distance  traveled  by  the  signal  from  the  platform  A  to  the  target 
plus  the  distance  from  the  target  to  point  B  (i.e.,  Rp  +  Rt).  Direc¬ 
tion  2  is  the  direct  distance  RA  from  the  platform  A  to  platform  B. 

The  signal  along  direction  1  is  received  at  Tx .  The  signal 
along  direction  2  is  received  at  T2 .  The  time  difference  is  AT=T2-TX. 
If  the  signal  is  a  response  from  the  target,  the  time  T2  includes  any 
fixed  and/or  random  time  delays  introduced  by  the  target  transponder. 
AT  has  an  associated  measurement  error,  dAT,  whose  average  value  is 
zero  and  whose  standard  deviation,  cdiT,  is  known. 

For  this  case  one  can  show  from  analytic  geometry  that  the 
target  location  will  lie  on  an  ellipse  defined  by  AT  and  R1  (see 
figure  8b).  The  ellipse  has  foci  at  (Ri-C0 AT/2,0)  and  (CoAT/2,0), 
and  vith  eccentricity  equal  to  -1+C0AT/Ri,  where  Co  is  the  speed  of 
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T  (x.y) 


A  (0,0)  B  (R|  ,0) 


Figure  8a.  Case  3:  Signal  Direction  of  the  Time  Measurements 
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light.  Geometrically,  target  location  will  be  the  intersection  of 
the  relative  bearing  line  n  +  0  -  (or  -  +t)  vith  the  ellipse. 
Analytically,  target  location  (x,y)  is  obtained  by  solving  the 
simultaneous  equations 


TanCe-^)  -  y/x  (  or  Tan(  ♦<-♦*)  -  -y/d^-x)  ) 
and 

RP  +  Rt  -  Ri  *  C»AT  • 

The  partial  derivatives  of  (x,y)  vith  respect  to  the  measurement 
variables  are  given  in  appendix  C.  The  Trace  and  Det  are  obtained  by 
substituting  these  derivatives  into  equation  (13).  Evaluation  gives 


Trace  « 


aBp'(ode  *<’d,*><1*cos®)  *  C>dCT  *  I l+cos(  ♦. -dj ) J 1  Ojr* 

(l+cos«)J 


and 


Det 


Rp<»d9*  "d,1)|Co<,dffl'  *  | 

(l+coso)J 


vhere  a  ■  -  ©• 

Figure  9  is  a  plot  of  CEP  contours  for  the  Case  3  geometry  vith 
-  0  deg,  -  ISO  nmi  and  vith  standard  deviations  ed9  -  5  deg, 
odt  *  0  deg,  odAT  -  10' 6  sec.  Platform  locations  A  and  B  are 

i 

marked  on  the  axes.  Similar  to  Case  1,  the  graph  is  to  be  interpreted 
as  follows:  Select  a  point  (x,y).  Linearly  extrapolate  (as  a  first 
approximation)  betveen  curves  to  estimate  the  curve  passing  through 
(x,y).  The  contour  number  associated  vith  that  curve  corresponds  to 
the  CEP  for  point  (x,y). 
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X  (Nautical  Miles) 

Figure  9.  Case  3j  CEP  Contours  for  One  Bearing  Angle 
and  a  TDOA  (Ellipse) 

a<te  -  5  deg,  aitT  -  10*6  sec,  o  -  0  deg,  R.  »  150  nmi 


As  discussed  earlier,  the  square  root  of  the  Trace  eay  be  used  to 
qualitatively  describe  the  CEP.  In  the  above  equation,  division  by 
zero  occurs  vhen  cos  (a)-  -1.  This  condition  occurs  along  the  hori¬ 
zontal  axis  of  figure  8  on  the  interval  between  A  and  B  (vhen  6*0, 

♦t  *  180,  or  0  *  180,  +t  *  0).  For  Case  3  geometry,  this  interior 
interval  on  the  horizontal  axis,  (and  a  small  band  around  it)  is 
called  a  "black  hole"  region.  Along  this  interval  the  CEP  is 
infinite.  Targets  near  this  interval  have  extremely  poor  location 
measurement.  In  contrast,  targets  outside  of  A  and  B  along  the 
horizontal  axis  are  able  to  be  located  with  good  accuracy. 

An  example  of  Case  3  application  is  identification,  friend  or 
foe  (IFF),  interrogation  of  targets  from  a  (location  A)  ground  or 
airborne  station. 
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SECTION  6 
CONCLUSIONS 


This  paper  describes  an  algorithm  for  calculating  Circular  Error 
Probability  (CEP).  CEP  relates  measurement  errors  to  the  variation 
of  a  calculation  that  is  computed  from  those  measurements.  As  such, 
it  is  a  useful  tool  for  assessing  the  accuracy  and  sensitivity  of  any 
calculated  quantity  to  errored  inputs.  CEP  is  frequently  used  with 
target  location  estimates. 

A  mathematical  foundation  is  provided  for  the  CEP  estimate,  and 
the  algorithm  is  shown  to  be  extremely  accurate.  It  differs  from  the 
true  CEP  estimate  by  less  than  one  percent  on  average  and  has  a 
maximum  error  of  1.5  percent. 

Three  common  cases  of  CEP  application  are  discussed.  All  deal 

with  assessing  the  accuracy  of  a  target  location  in  a  bistatic 

environment.  Examples  of  CEP  contours  are  generated  for  the  special 

conditions  Ri  =  150  nmi,  ^  *  0  deg,  odR  =  0  nm,  od+  =  0  deg. 

i  i 

Regions  of  poor  target  detectability  are  indicated. 
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APPENDIX  A 

DERIVATION  OF  THE  SIMPLIFIED  CEP  FORMULA 


CEP  is  defined  as  that  value  of  R  that  satisfies  the  integral 
equation 


p(u,v)dudv  =  0.5 


u2+v2<  P. 


2 


(Al) 


where  p(u,v)  is  the  joint  Gaussian  distribution  of  variables  (u,v) 
with  mean  (0,0)  and  variances  (ou2,ov2)  » 


p(u,v) 


2  n®u  °v 


Htf- tttn 


(a2) 


Ve  make  the  following  change  of  variables: 


Let 

aL-  Maximum!  ou,<xv]  , 

os*  Minimum!  «Tu,crv]  , 

(A3) 

w  =  ffs/oL,  0<  w  <  1 

1 

and  define 

s  *  u/aL  , 

(A4) 

t  -  v/aL  . 

If,  for  the  sake  of  discussion, 

we  suppose  that  oL  * 

<ru  (the  same 

argument  will  hold  for  oL-ev), 

then  integral  (Al)  can 

be  reexpressed 

as 
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Symbolically,  we  write  (Al)  in  the  form 

Fpi, 

L  ffL 

to  emphasize  the  fact  that  integration  of  (A5)  results  in  an  implicit 
relation  between  (CEP/ a.  )2  and  w. 

In  general,  it  is  not  always  possible  to  obtain  an  unambiguous 
solution  from  an  implicit  relation.  (Remember  that  our  object  is  to 
solve  (A6)  for  CEP).  However,  if  it  can  be  shown  that  (CEP2/aL2)  is 
an  increasing  function  of  w,  then  the  Inverse  Mapping  Theorem  of 
Calculus  argues  that  we  may  uniquely  solve  equation  (A6)  for 
(CEP2/oL2)  as  a  function  of  w. 

We  show  that  (CEP2/aL2)  is  an  increasing  function  of  w  by 
proving  that  its  derivative  with  respect  to  w  is  positive. 
Differentiating  (A6)  we  obtain 

3F  3F  d(CEP2/o2) 

3w  +  3TCEP2/<72)'  dw  =  0  .  (A7) 

Solving  for  the  derivative  in  (a7)  gives 


d(CEP2/o2)  3F/3v 

dw  3F/3(CEP2/<£) 


(A8) 
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The  positiveness  or  negativeness  of  the  partial  derivatives  on  the 
right-hand  side  of  (A8)  can  be  obtained  by  the  following  arguments: 

1.  As  (CEP2 /or2)  increases,  the  circle  area  over  which  integration 
takes  place  increases.  Therefore,  the  value  of  the  integral  F 
will  also  increase  for  each  fixed  v.  Ve  conclude  that 

3F 

3TCEP2/ff.2)  >  0  ' 

2.  As  w  increases,  the  limits  of  integration  remain  constant  but 
the  integrand  decreases.  Therefore,  the  value  of  the  integral 
(i.e.,  F)  will  decrease  for  each  fixed  (CEP2 /a2). 

p 

Ve  conclude  that 


The  above  inequalities  guarantee  that  the  derivative  in  (A8)  is 
positive.  Thus  (CEP2/o2)  is  a  monotonically  increasing  function  of 
w. 

The  Inverse  Happing  Theorem  of  Calculus  can  now  be  applied.  It 
states  that  ve  can  (in  theory)  uniquely  solve  equation  (A6)  for 
(CEP2/c  2)  as  a  function  of  v.  Thus  for  some  function  h(v)  we  will 
have 

h(w)  and  CEP  -  crLJh(wy  .  (A9) 

Equation  (A9)  gives  an  explicit  representation  for  the  form  of  the 
CEP  solution  as  a  function  of  w. 
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The  function  4h(v)  can  be  tabulated  numerically  from  equation 
(Al)  by  selecting  a  value  for  v  and  then  solving  for  the  value  of 
CEP/ol  which  will  make  the  integral  equal  0.5.  Figure  5  shows  the 
plot  of  the  tabulated  point  pairs  (w,CEP/oL).  One  notes  that  on  the 
interval  0  <  w  <  0.5  the  curve  follows  a  quadratic  pattern,  while  on 
the  interval  0.5  <  w  <1  the  curve  is  essentially  linear. 


The  simplified  approximation  formula  for  CEP  is  nothing  more 
than  a  reflection  of  this  observation.  It  involves  a  quadratic  fit 

to  the  curve  for  ^lh(w)  on  the  interval  0  <  w  <  0.5  and  a  linear  fit 

h(w)  on  the  interval  0.5  <  w  <1,  namely, 


^lh(w)  = 


(  0.67  ♦  0.8w2 
l  0.59(1  +  w) 


if  0  <  w  <  0.5 
if  0.5  <  w  <  1 


(A10) 


Equation  (A10)  is  an  excellent  approximation  to  CEP/aL.  On 
average,  the  error  to  the  true  curve  is  less  than  IX  over  the 
interval  0  <  w  <1. 


APPENDIX  B 

THE  COVARIANCE  OF  (dx,dy) 


The  covariance  of  a  vector  variable  v  with  mean  0  is  defined  by 

Cov(v)  *  E(vv*)  ,  (Bl) 

where  E(*)  is  the  expected  value  operator  and  v*  means  the  transpose 
of  v.  The  operator  E  is  linear  and  has  the  property  for  constants 
a,b  and  variables  vx,v2  that 

E(av1  +  bv2)  =  aE(vx)  +  bE(v2)  .  (B2) 

If  v  is  of  dimension  N,  the  covariance  of  v  will  be  an  NxN  matrix. 

Let  v  -  \.dx,dy)*  and  q  =  (dmx  ,dm2  .dm3  )* .  From  equation  (5)  of 
section  2  we  are  given  the  relation 

v  =  Mq  ,  (B3) 

where  M  is  a  matrix  of  constants. 

The  following  :<-;.j,u  tat  ions  now  give  us  the  desired  relationship: 

vv*  ■  (Mq)(Mq)*  *  Mqq*M*,  (from  property  of  transpose) 

E(vv*)  *  E(Mqq*M*'.  ME(qq*)N*,  (from  B2) 

Cov(v)  *  ME(qq*/  -•  MCov(q)M*.  (from  Bl,  def.  of  covariance) 


APPENDIX  C 

CEP  DERIVATIVE  COMPUTATIONS  FOR  CASES  1,  2,  3 


Case  1.  *1^  »+t  »9)  are  measured  with  error  variances  (od+,  <xdR  » 

i 

°df  *  CTde ) •  Geometrically,  target  location  (x,y)  is  obtained  as  the 

t 

intersection  of  two  bearing  lines.  From  the  triangle  relationships  in 
the  geometry  of  figure  3,  the  equations  from  which  target  location 
can  be  obtained  analytically  are 

Tan(&-*A)  =  y/x  and  Tan(*t-fA)  =  -y/(RA-x)  . 

The  partial  derivatives  of  x  and  y  are  given  by 


ax 

Rt  cos  (  0-  ^  )  -Rp  cos  ( *t  -  ) 

3x 

cos(9-^1)sin(+t-*i) 

sin(a)  ' 

sin(+t-0) 

3x 

Rtcos(0-^) 

3x 

-RjSMft-^) 

“  f 

— 

—  1 

sin(4t-9) 

36 

sin(ft-6) 

3y 

Rtsin(6-4i  )-Rpsin(#t-4i ) 

9y 

cosC©-^ )sin( ) 

3*i 

"  t 

sin(a) 

dRi 

sin(  -9) 

3y 

F^cosCS-^ ) 

3y 

-Rtsin(0-^) 

3*t 

sin(*t-0) 

30 

sin(*t-0) 

where  a  » 

*  *t  -  e. 
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Case  2. 


(  0,  AT ,  Rj ,  ^  ,  4t )  are  measured  with  error  variances 
( ade » ffdiT '  °<ir  *ffd*  )•  The  tRr8et  location  will  lie  on  an 

i  i  t 

hyperbola  defined  by  AT  and  Ri .  The  hyperbola  has  foci  at 
(-CoAT/2,0)  and  (Ri+Co AT/2,0),  and  eccentricity  equal  to  Ri/C0AT, 
where  CQ  is  the  speed  of  light.  Geometrically,  target  location  will 
be  the  intersection  of  the  relative  bearing  line  n  +  0  -  ^  (or 
~  )  yith  the  hyperbola.  The  equations  from  which  the  target 

location  (x,y)  can  be  obtained  analytically  are 


TantO-^)  -  y/x  (  or  Tan(*t-^)  «  -y/(Rrx)  ) 

and 

Rp  +  ^  -  Rt  =  CoAl  . 

The  partial  derivatives  of  x  and  y  are  given  by 


3x  Rp[sin(*t-+i)-sin(0-^.  )] 

l-cos(  a) 

3x  RpIsinC^-^l-sinCe-^)) 

30  1-cos  (a) 

3y  RplcosC^-^i-costO-^)! 

3^  ~  1-cos  (a) 

3y  RplcosC^-^i-cost  0-^)1 

sx  —  '  - 

30  l-cos(a) 


3x 

w, 

3x 
3  AT 


C0  008(0-^) 
1-cos  (a) 


3y  fl+cos(* 

3R,  “  ”  L  l-cos(a)  ~Jsin(e'^)  * 


3y  Cosin(0-+i) 

s  ■  ■ 

3  AT  l-cos(a) 


where  a  -  -  0  and  Co  is  the  speed  of  light. 


46 


Case  3.  ($t ,  AT,Ra , ^ )  are  measured  with  error  variances  (ed+  , 

t 

°«jAT»adR  >  °d*  )•  The  target  location  will  lie  on  an  ellipse 

i  7  i 

defined  by  AT  and  R^  The  ellipse  has  foci  at  (Ri-C0 AT/2,0)  and 
(Co AT/2,0),  and  with  eccentricity  equal  to  -1+C0AT/Ri,  where  Co  is 
the  speed  of  light.  Ge metrically,  target  location  will  be  the 
intersection  of  the  relative  bearing  line  I!  +  0  -  (or  -  ♦t ) 
with  the  ellipse.  The  equations  from  which  target  location  (x,y)  can 
be  obtained  analytically  are 

Tan(e-^)  «  y/x  (  or  Tan(*t-^)  =  -y/fR^x)  ) 
and 

R  +  R  -  R.  »  C  AT  . 

p  t  i  o 

The  partial  derivatives  of  x  and  y  are  given  by 


ax 

%lsin(#t-#1)+sin<e-#i)l 

dx 

ri+cos(*  -♦  )] 

a?x  - 

l+cos(a)  ’ 

3 

L  l+cos(a)  Jcos<*-*i>  » 

ax 

Rplsin(+t-4i)+sin(0-4i )] 

cocos(e-41) 

ae 

> 

l+cos(o) 

3 

a  at 

s  # 
l+cos(a) 

ay 

Rp  [  cos(  )+cos(  )  i 

ay 

p+cos(#t-4i)1 

l+cos(a)  ’ 

**  3 

■  [  l+cos  ( a)  Jsi”(6-*i)  . 

ay 

Rp(cos(^t-4i)+cos(e-*1)J 

ay 

C^since-^) 

ae 

t 

l+cos(a) 

a 

a  at 

g  f 

l+cos(a) 

where 

a  ■  ♦t-®  and  Co  is  the  speed  of  light. 

APPENDIX  D 

TARGET  LOCATION  WHEN  DISTANCE  BETWEEN  PLATFORMS  IS  NOT  KNOWN 


An  interesting  but  more  complex  case  occurs  if  one  assumes  that 
the  distance  Rt  between  the  platforms  is  an  unknown.  In  this  case, 
two  bearing  angle  measurements  and  a  TDOA  of  Case  3  type  (ellipse) 
will  provide  a  solution  not  only  for  target  location,  but  also  for  Rt 

Relationships  between  the  sides  and  angles  of  triangle  A,B,T  in 
the  geometry  of  figure  2  are  given  by 

R,  R>  Ri 

-  a  -  _  -  (The  Law  of  Sines) 

sinC^-^)  sin(0-*.)  sin(a) 

y  «  -(Ri-x)tan(  ♦t-f1 )  and  y  =  xtante-^)  , 

A. 

Rp  +  Rt  -  Ri  ’  CoAT  ‘ 

The  solution  of  these  simultaneous  equations  for  x,y  and  Ri  is 

Co  ATsin(o) 

Ri  “ 

sinCS-^ )  +  sin(#t-ft)  +  sin(a) 

Co  ATcos  (  e- ^ )  s  i  n  ( ) 

x  *  - 

sin(0-^i )  +  sin(4t-^)  +  sin(o) 

CoATsin(0-4i)sin(^t-#i) 

y  =  - 

sinCe-^)  +  sin(4t-f1 )  +  sin(a) 

where  a  ■  -  0  and  Ce  is  the  speed  of  light. 
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The  partial  derivatives  of  x,  y  and  R  are  given  by 

8x  Rtco8(S-4i)IUcos(#t-^i)l 

3?t  *  D  * 

dx  Rp  [  l+sin(  )sin(^-4i  )+cos(  4t-4t )  1 

—  » - _______  f 

36  D 

3x  Coc°s(^-*t)8inC»t-»i) 

3  AT  D 

3y  Rtsin(6-4i)[l+cos(^-4i)] 

8?t  *  D  ’ 

3y  RpSin(+t-4.  )ll-cos(6-ti)] 

—  »  -  * — •-  . "  _  .  , 

36  D 

3y  Cosin(6-ii)sin(ft-4i) 

3  AT  0 


*1 

36 

3  AT 


Rt  (sin(a)cos( ♦t-*i)+sin(  ^-^J+sinf  6-^)1 

oiTKTvv 

Rp{sin(a)sin(6-4i)+cos(6-4i)+cos(#t-*i)J 

Dcos(9-^) 


C0sin(a) 


where  D  -  sin(+t-tA)  +  sin(6-^)  -  sin(a)  , 

a  ■  4t  '6  and  C0  is  the  speed  of  light. 
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